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ABSTRACT 



The spatial distribution of satellite galaxies around host galaxies can illuminate the relationship between satellites and dark matter 
subhalos and aid in developing and testing galaxy formation models. Previous efforts to constrain the distribution attempted to elimi- 
nate interlopers from the measured projected number density of satellites and found that the distribution is generally consistent with 
the expected dark matter halo profile of the parent hosts, with a best-fit power-law slope of -1.7 between projected separations of 
~ 30/;"' kpc and 0.5/i"' Mpc. Here, I use the projected cross-correlation of bright and faint galaxies to analyze contributions from 
satellites and interlopers together, using a halo occupation distribution (HOD) analytic model for galaxy clustering. This approach 
is tested on mock catalogs constructed from simulations. 1 find that analysis of Sloan Digital Sky Survey (SDSS) data gives results 
generally consistent with interloper subtraction methods between the projected separations of lOIr' kpc and 6.3/;"' Mpc, although 
the errors on the parameters that constrain the radial profile are large, and larger samples of data are required. 

Key words, cosmology: theory - dark matter - galaxies: formation ~ galaxies: structure - galaxies: fundamental parameters 



1. Introduction 

In the hierarchical assembly of dark matter (DM) halos, progen- 
itor halos merge to form larger systems. Some DM halos survive 
accretion and form a populati on of subhalos within the virial 
radius of the larger halo (e.g.. iGhigna et al.]| 1 998t iKIvpin et alj 
. 119991) . Baryonic material cools and forms stars in the center 
' of some halos and subhalos, resulting in galaxies and satel- 
lite galaxies. Satellite galaxies are biased tracers of the poten- 
tial well of their host DM halos, and satellite dynamics have 
' been used to provide constraints on the total mass distribution in 
' galactic halos (e.g ., Zaritskv & White 1994; Zaritskv et al. 1997; 
'Prada et al."2003'; van den Bosch et alj2004j;iConrov et al.i2005t: 
iFaltenbacher & Diemand 2006) and to test the nature of grav- 
, ity (Klvpin & Prada 20071). In addition, the spatial distribution 
of the satellite galaxy population reflects the evolution of dwarf 
galaxies and the mass accretion history of the halo. For exam- 
ple, observations of satellites around early-type galaxies suggest 
that the distr ibution lie s along the m ajor axis of the light distri- 
butio n (e.g.. ISales & L ambas 200l iBrainerdlllOOl lYang et al.1 
120061; lAzzaro et al.fl20b 7). while in DM simulations the angu- 
lar distribution of subhalos follows the shape of the DM halo, 
which is indicative of infall of subhalos along filaments (e.g., 
IZentner et ai]|2005l; iLibeskind et al.ll2005l) . 

The radial spatial bias of satellite galaxies is another impor- 
tant indicator, providing a simple observable for testing galaxy 
formation models. Simulations of galaxy cluster-sized halos sug- 
gest that the distribution of DM subhalos has a radial pro- 
file that is less concentrated than the DM density distribution 
of the host halo at small radii (within ~ 20 - 50% of the 
virial radius) but has a profile similar to the DM halo at larger 
radii (Ghignaetal. 1998; CoKn et al. 1999; Ghignaetal. 2000; 
Springel et al. 2001; De Lucia et al. 2004; Diemand et al. 2004; 
Gao et al.M2004; .Nagai & Kravtsov, ,2005; .Maccio et al.,,2006i) . 



Simulations that include baryons, star formation and cooling, 
however, show that the distribution of galaxies associated with 
subhalos has a steeper inner profile than the subhalo distribu- 
tion, both at cluster and at galaxy scales, and these distributions 
are in general agree ment with observations (iNagai & Kravtsovl 
l2005t iMaccio et al ^'2006). For samples of subhalos selected by 
dark matter mass, objects found near the halo center have lost a 
greater percentage of their dark matter mass due to tidal stripping 
than objects near the virial radius. As result, the concentration 
of the distribution of subhalos is lowered. Stellar mass-selected 
samples of satellite galaxies are resistant to this elfect since bary- 
onic components are located in the centers of dark matter subha- 
los and are tightly bound. 

Observations of the Local Group dwarf population suggest a 
distribution more radially concentrated tha n that of DM subhalos 
in dissipationless numerical simttlation s jKravtsov et alj 120041; 
iTavlor et al.ll2004t IWillman et al.1120041 ). At large radii, this dis- 
tribution would also be inconsistent with the DM halo profile. 
This strong radial bias, however, may not apply to brighter satel- 
lites - such as the Magellanic Clouds - and may be limited to 
the faintest dwarfs, related to the physics of the formation of 
the smallest dwarf galaxies jKravtsov et alj2004tlDiemand et al.l 
2005) or incompleteness in observing faint objects at large radii. 

Constraints on the radial distribution of galactic satellites be- 
yond the Local Group have been limited by the effect of inter- 
lopers - projected objects contaminating the sample of satellite 
galaxies. Previous attempts to constrain the satellite distribution 
have used galaxy redshift surveys and various methods of sub- 
tracting interlo pers from samples of faint galaxies near blight, 
~ Lt, galaxies. Ivan den Bosch et aTl (12005! ) used data from the 
Two Degree Field Galaxy Redshift Survey (2dFGRS), but con- 
cluded that the data does not allow meaningful constraints on 
the radial distribution at small radii due to the incompleteness of 
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close pairs in the survey, although the data is consistent with the 
dark matter distribution at large radii. In an independent analysis, 
[Sales & Lambas (2005) fit a power-law slope to the satellite dis- 
tribution for projected radii between 20 and 500 kpc, finding 
a slope shallower than expected for the DM halo with a signif- 
icant dependence on morphological type of the pare nt galaxies 
but using no interloper subtraction. iChen et alj ( 12006 ) used data 
from the Sloan Digital Sky Survey (SDSS), testing methods of 
interloper subtraction on numerical simulations to show that the 
radial profile of satellites of isolated ~ L, galaxies is steeper 
than the subhalo profile in simulations and possibly consistent 
with being as steep as the DM profile. 

An alternate approach to this problem is to use the projected 
cross-correlation of bright and faint galaxies. The correlation 
function is the statistical measure of the excess probability over 
a random distribution of finding pairs of objects at separation, r. 
It can be calculated analytically from the mass function of halos, 
the spatial distribution of galaxies within halos, the halo bias, 
and the probability of finding a number of galaxies in a halo 
of mass, M, parameterized by the halo occupation distribution 
(HOD ; Ma&Frv 2000; Seljak 2000; Scoccimarro et al. 2001; 
iBerlind & We inberg 2002). Many recent works have used such 
analyses of autocorrelated samples of galaxies in the SDSS in 
order to constrain the relation between distributions of galax- 
ies and dark matter and to constrain cosmological parame- 
ters (e.g., Zehavi et al. 2002, 2004, 2005; Abazaiian et al. 2005; 
IZheng et al.ll2007l:lTinker et al.li2008n . 

A halo model analysis of the clustering of galaxies can con- 
strain both the satellite and interloper populations, without the 
need to cut the bright galaxy sample by environment or use sub- 
traction methods in order to eliminate interlopers. Here, I use the 
cross-correlation of bright and faint galaxies, as well as the au- 
tocorrelation functions of bright and faint galaxies, to constrain 
the spatial distribution of faint satellites around bright galaxies, 
using data from the SDSS spectroscopic sample. The ability of 
the data to constrain the spatial distribution of satellites with this 
method is tested with mock data sets from simulations populated 
with galaxies by a realistic HOD. Best-fit parameters and errors 
are characterized using Monte Carlo Markov chains (MCMC). 
Despite the degeneracies between parameters, constraints on the 
spatial distribution appear without systematic bias. Errors on the 
parameter estimates, however, are large in the current data set, 
and larger datasets are required. 

The paper is organized as follows. The SDSS data used in the 
analysis are discussed in Section |2l while the coiTelation func- 
tion and associated statistical and systematic eiTors are discussed 
in Section [3] The HOD analytic model is detailed in Section 
m and modeling with MCMC and populated simulations is dis- 
cussed in Section|5] The main results are presented in Section|6l 
Conclusions are summarized in Section|2l Throughout, I assume 
flat ACDM cosmology with D.^ = 0.3, Qb = 0.04, erg = 0.9, and 
h = 0.7. 

2. The Sloan Digital Sky Survey 

The Sloan Digital Sky Survey (SDSS) will image up to 10"^ 
deg^ of the northern Galactic cap in five bands, u,g,r,i,z, 
down to r ~ 22.5 and acquire spectra for galaxies with r- 
band Petrosian magnitudes r < 17.77 and r-band Petrosian 
half-light surface brightnesses //so < 24.5 mag arcsec"^, us- 
ing a dedicated 2.5m telescope at Apache Point Observatory in 
New Mexico (Fukugita et al. 1996; Gunn et al. 1998; York et al.' 
2000; Hoggetal. 2001; Smith et al. 2002; Strauss et al. 2002; 
Blanton et al., ,2003bi: ,Gunn et ah i2006i; iTucker et ah i2006) . 



An automated pipeli ne measures the redsh i fts and classifies 
the reduced spectra (Stou ghton et all 120021; iPier et all 120031; 
Ivezic et al. 2004, D. J. Schlegel et al. 2009, in preparation). 

For this analysis, I use a subset of the spectroscopic main 
galaxy catalog available as Data Release Five (DR5) and in- 
cluding all of the data avai lable as of Data Release Four (DR4) 
dAdelman-McCarthv et al. 2006). I refer to this sample hereafter 
as DR4-I-. Because the SDSS spectroscopy is taken through cir- 
cular plates with a finite number of fibers of finite angular size, 
the spectroscopic completeness varies across the survey area. 
The resulting spectroscopic mask is represented by a com bina- 
tion of disks and spherical polygons (Tegm ark et al.l l2004). Each 
polygon also contains the completeness, a number between and 
1 based on the fraction of targeted galaxies in that region which 
were observed. I apply this mask to the spectroscopy and in- 
clude only galaxies from regions where completeness is at least 
90%, for a completeness-weighted area of 5104 deg^. I use r- 
band magnitude s in DR4-I-, built fro m the NYU Value-Added 
Galaxy Catalog dBlanton et al.ll200"5h . normalized to h-l, such 
that 

Mr = Mo.i, - 51ogio/z, (1) 

where Mq.i, is the absolute magnit ude K-corrected to z=0.1 
(kcorrect v3 . 4) as described in Blanton et al.1 (l2003ah . 

A volume limited sample of galaxies down to Mr = -18 con- 
tains galaxies out to redshift z - 0.048, with a median galaxy 
redshift of z=0.038. The redshift limit is a trade-off" with the 
magnitude limit. I create a 'bright' sample and a 'faint' sam- 
ple, defined such that bright galaxies have r-band magnitudes 
between -21 and -20 and faint galaxies have magnitudes be- 
tween -19 and -18. The bright galaxy sample contains galaxies 
that are near L*, while the faint galaxy sample contains galax- 
ies that have luminosities similar to bright satellites such as the 
Magellanic Clouds. By excluding objects between Mr - -19 
and Mr - -20, the samples are fairly analogous to the samples 
used in previous works t hat cons trained the spatial distribution 
of satellite galaxies (e.g.. lChen et al. 2006). The corresponding 
number density of objects for these samples is calculated from 
the data and shown in Table [1] along with the total number of 
objects in each sample. 

In addition to the faint and bright samples, a luminosity- 
threshold sample of Mr < -21 and a luminosity-threshold sam- 
ple of Mr < -19 are created from the same volume-limited sam- 
ple. These samples are necessary in order to constrain the masses 
of the halos in which bright and faint galaxies are found, as fur- 
ther described in Section |4] 

Table 1. Number Densities of Volume-Limited Samples of 
SDSS Galaxies 



My 




N 


-21 to -20 ('bright') 


0.00392 


5863 


-19 to -18 ('faint') 


0.01137 


16992 


< -19 


0.012549 


19085 


< -21 


0.00056 


846 



3. Estimating the Correlation Function 

The correlation function measures the excess probability of find- 
ing a pair of objects at some separation, comparing counts of 
pairs of objects in catalogs of real objects (DD) and in catalogs 



Jacqueline Chen: The Galaxy Cross-Correlation Function as a Probe of the Spatial Distribution of Galactic Satellites 



3 



of objects with randomized positions (RR) and pairs where one 
object is in the real catalog and its counterpart in the random 
catalog (DR or RD)[3The total number of pairs is normalized by 
the total number of objects in a sample, such that 



DD = 



number of pairs 



(2) 



NdNd 

Using the iLandv & Szalavl (Il993h estimator, the cross 
correlation as a function of projected separation, r^, and line 
of-sight separation, n, is 



DD-DR-RD + RR 



RR 



while the autocorrelation is 
DD - 2DR + RR 



RR 



(3) 



(4) 



The parallel and perpendicular components of the pair sepa- 
rations are distinguished in the data by 



si 



and 



2 2 



(5) 



(6) 



where s = Vi - V2, 1 = ^(vi + V2) and Vi and V2 are the redshift 
positions of the two objects in a pair (Fisher et al. 1994). 

From the correlation function, the prq/ecfeii correlation func- 
tion is 



Wpirp) 



^{rp,7T)dn. 



(7) 



where n is integrated to tt = 40/i"' Mpc, large enough to in- 
clude all significant clustering signal. The projected correlation 
is measured in logarithmic bins for projected radii between 10 
/i"' kpc to 6.3 /i"' Mpc, with 14 total data points. 

Statistical eiTors in the correlation are estimated via jack- 
knife sampling (see, e.g., Lupton 1993). In the SDSS data, the 
survey area is divided into 205 equal area samples using the hi- 
erarchical pixel scheme SDSSPijfl which represents well the 
rectangular geometry of the SDSS stripes. The covariance matri- 
ces are estimated by the covariance between jackknife samples 
- samples where one of the 205 equal area subsamples is omit- 
ted. In order to appropriately estimate the covariance matrix, the 
number of jackknife samples must be significantly greater than 
the number of bins in which the correlation function is estimated 
dHartlap et al.l 1200 7). Typically, each jackknife sample is con- 
structed to be larger than the largest separation measured in the 
corresponding coiTelation function, and the number of jackknife 
samples is at least equal to the square of the number of bins. 
With 14 bins, 196 is the minimum number of jackknifes sam- 
ples. The data sample uses slightly more, 205. The typical jack- 
knife region, then, has a comoving volume of ~ 19"'/z^Mpc"^. 
The projected cross-coiTelation and corresponding autocorrela- 
tion functions are shown in Fig.[T] 



' For the SDSS samples, the random catalogs contain at least twenty 
times the number of objects in the real catalogs. In order to create a 
random catalog from a data catalog, the objects in the data catalog are 
sampled with replacement - such that each data point may be sampled 
several times - and then assigned a randomized position on the sky. 

^ See,http://lahn[iu.phyast.pitt.edu/~scranton/SDSSPiX| 
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Fig. 1. The measured projected correlation functions in the 
SDSS data, where the cross-correlation of bright and faint galax- 
ies is indicated by filled circles, the bright autocorrelation by 
open circles, and the faint autocorrelation by filled triangles. 
Bright galaxies have r-band magnitudes between -21 and -20, 
while faint galaxies have r-band magnitudes between -19 and 
-18. 



In the SDSS spectroscopic survey, no single pointing of the 
telescope can measure the spectra of objects that are separated 
by less than 55", the fiber collision distance. At the median red- 
shift of the sample, 55" is 30 /j"' kpc, well above the minimum 
separation of 10 h^^ kpc in the data. To achieve the minimum 
separation, I employ a catalog corrected for collisions. For pairs 
of objects separated by less than the fiber collision distance in 
the target selection catalog, the unobserved target is assigned the 
redshift of the observed object in the pair. This could bias the 
sample by overestimating pairs, but significant error is unlikely 
given that pairs at those separations are extremely correlated (as 
in Fig.O. 

Previous studies have described systematic errors with angu- 
lar depend ence (Masjedi et al. 2006; Mandelbaum et al. 2006^ 
lAdelman"-McCarthv et alj 120081) . Issues with sky subtraction 
may lead to underestimation of pairs by ~ 5% with angular sep- 
arations between 40 - 90" and an overestimation of pairs at an- 
gular separations less than 20 - 30". This effect is most acute for 
pairs of objects that consist of an apparently bright object and an 
apparently faint object (Rachel Mandelbaum 2006, private com- 
munication). 

In the cross-correlated sample, the estimated systematic er- 
ror of 5% for pairs with angular separations between 40" and 
90" is smaller than the statistical eiTors in the coiTelation func- 
tion, ~ 10%, indicating that such errors do not affect the results. 
For the unconstrained overestimate of pairs with angular sepa- 
rations less than 20", only the first bin of the estimated cross- 
coiTelation (containing objects with projected physical separa- 
tion of less than 16 kpc) is affected. In this bin, 23% of the 
pairs could be possibly affected, while the fractional statistical 
error in the correlation function is 15%. Given that only one of 
the data points is affected, it seems unlikely that the parameter 
estimates are biased by this systematic. However, I test the pos- 
sible effects of this error in Sections [5.51 and [63] 
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4. Halo Model of Correlation Functions 

The correlation function can be modeled using the halo model - 
assuming that all galaxies and dark matter particles live in dark 
matter halos. First, the correlation function is decomposed into 
two parts: a one-halo term and a two-halo term. 



^(r) = fih(r)+^2hW, 



where the terms correspond to contributions from galaxy pairs 
found in the same dark matter halo and galaxy pairs found in 
different halos. At small separations, r, the one-halo term dom- 
inates the coiTelation function, while at large separations, the 
two-halo term dominates. 

All galaxies are identified as either a central galaxy, found 
at the center of a dark matter hal o, or a s atellite galaxy, found 
elsewhere in the halo. Following [Tinker e t al. (2005, see their 
Appendix A & B), the one-halo term is decomposed into two ad- 
ditional terms: (1) a central-satellite term where one galaxy of a 
pair is the central galaxy and the other is a satellite of the central 
galaxy and (2) a satellite-satellite term where both galaxies in 
a pair are satellite galaxies in a halo. As in Berlind & Weinberg! 
(120021 Eqn. 11), the one-halo term can be written. 



1 



1112 



:J-dM^{N,,,,)Mj^F'{^), (9) 



where hi and h2 are the number densities of galaxies, and here 
and throughout this section ' 1 ' and '2' are used to denote generic 
samples. dn/dM is the mass function of halos, and (A^pair)M is the 
average number of pairs of objects in a halo of mass, M. 7?vh is 
the virial radius of a halo of mass, M, and F(x) is the average 
fraction of galaxy pairs in such a halo with separation less than 
r, and F'(x) is the derivative with respect to r. The virial radius 
of a halo is defined as the radius at which the mean mass density 
enclosed is 200 times the mean matter density of the universe. 

To separate the central-satellite term from the satellite- 
satellite term, for the cross-coiTelation function. 



{Np^u)mF' (x) = ((A^cen,bA^sat,f) + (A^cen.f A^sat,b))^'cen-sat(-^) 

and for the autocoiTelation, 



(10) 



{Np,,,}mF' (x) 



(11) 



where centrals are denoted by 'cen,' satellites are denoted by 
'sat' and bright and faint objects are labeled 'b' and 'f.' f'^en-sat 
is proportional to the volume-weighted halo density profile and 
^sat-sat is proportional to the convolution of the halo density pro- 
file with itself. Den sity profiles are modeled with NFW profiles 
dNavarro et al.lll997|) arid the co nvolution is calculated analyti- 
cally as derived by fshethet all (2001) and Zheng (2004). The 
concentration of the dark matter halos is calculated following 
iBullock etal](l200lh from a simple expression for concentration 
as a function of the redshift of collapse for a halo. The galaxy 
distribution within the dark matter halos is parameterized by a 
scaling factor, such that C gai - /cdm- Finally, the mass function 
of halos from I Jenkins et aLl (12001 ) is employed. This fitting for- 
mula approximates the mass function of halos with an accuracy 
better than x 20% down to halos of lO'^/i^'M©, which is well 
below the range of masses used in the data and populated sim- 
ulations (see Sections |5] and |6]l. As mentioned previously, a flat 
ACDM cosmology with = 0.3, Qb = 0.04, as = 0.9, and 
h - 0.7 is used. 



The two-halo term is more easily calculated in Fourier space, 
where the power spectrum P(k) is the Fourier transform of the 
correlation function ^(r): 



(8) P,„{k) 



-r 



dM, 



dn 
'dM\ 



■{N)MMMur)y„{k,Mi) 



"2 J 



dMy 



dn 

-dm 



{N}MMM2,r)yg(k,M2) 



(12) 



where yg{k, M) is the Fourier transform of the halo density pro- 
file and {N)m is the average number of objects in a halo of mass, 
M. Common simplifications to this calculation have employed 
the large-scale bias, b{M), and approximated the matter power 
spectrum, P,„ik), with the linear power spectrum, Punik). Here, 
I use the dark matter power spectrum of iSmith et al.i (i2003i) and 
the scale-dependent halo bias such that 



[I + 1.17^„,(r)] 
[l+0.69^„,(r)] 



1.49 



2.09' 



(13) 



where the scale-independent bias is provided by [Tinker et alj 
see their Appendix A). They calibrate the analytic model 
of Sheth & Tormen ( 1999) with numerical simulations. In addi- 
tion, the assumption that the separation between halos is much 
greater than the scale of the halos leads to an overestimate in 
power at intermediate separations where the one-halo and two- 
halo terms are comparable - counting pairs in overlapping halos 
that would rightly be counted a s single halos . I use the spherical 
exclusion approximation from [Tinker et al.l ( l2005b . They only 
specify the two-halo term for the autocorrelation. For the cross- 
coiTelation, the geometric mean of the the two autocorrelated 
samples is used as an approximation, which performed well in 
tests on populated simulations. 

While the terms in the calculation for the halo mass func- 
tion, halo bias, and halo density profile are fixed, the number 
of galaxies and their spatial distribution in halos is parameter- 
ized and fit to the data. I use the HOD prescription for charac- 
terizing the galaxy bias, which has been s hown in simulations 
and semi-analytic models to be reasonable dBerlind et al.ll2003l : 
IZheng et"ani2005h . For luminosity-threshold samples (samples 
that include all objects brighter than a threshold luminosity) the 
average number of central galaxies per halo is a step function 
at a minimum mass, (A^cen) = 1 for M > and (A^cen) = 0, 
otherwise. For luminosity-binned samples (including all objects 
within a minimum and maximum luminosity), (A^cen) is a step 
function, with both a minimum mass and a maximum mass: 

<A^cen> = 1 forMmin <M < Mn,ax(cen)0 

For the average number of satellite g alaxies, a simple fu nc- 
tional form with three parameters is used dTinker et alj|2005b : 



(A^sat) 



^fexpf-^). 
Ml ^\ Ml 



(14) 



Using luminosity-threshold samples, simulations show that 
(A^sat) should "roll-off"' faster than a power-law for halo 
masses near the minimum mass, and HOD parameterizations 
that account for "roll-off" provide good fits to observations 



' iTinker et alj i2007l) suggest a more physical model for (A'cen) that 
takes into account the scatter between mass and luminosity. Tests com- 
paring the simple model for (A'cen) used in this work to their model, 
however, find no significant effect on the resulting correlation function. 
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dZehaviet al.ll2005h . In addition. IConrov et all (l2006h suggests 
that by adding a parameter for "roll-off," best-fits for the slope a 
may be driven to ff = 1 . However, for luminosity -binned samples 
- as the cross-correlation uses - parameterization of "roll-off" is 
not necessarily the same as is expected in luminosity-threshold 
samples. 

For each galaxy sample, there are three free HOD parameters 
~ Ml, a, and Mcut - and two fixed HOD parameters - Mmi,, and 
A^max(cen)- Afmin is Set by the observed number density of objects 
in the faint or bright sample such that 



«gal(f/b) - 



f 

Jm„ 



dn 

«A^cen (f/b)> + (A^sat (f/b)))— t/M 



dn 
dM 



dM 



dn 

{N,.,ifM)—dM 



(15) 



where (f/b) is used to note that the equation must be solved 
twice, once for the faint sample and once for the bright sample. 
Mniax(cen), however, cannot be set using those samples. Instead, 
it must be set by the number density of objects in samples that 
contains objects brighter than the objects in the corresponding 
bright or faint sample such that 



n 



gal (>f/>b) 



•J Ma 



dn 

((A^cen (>f/>b))+(A^sat (>f/>b))) — t/M, (16) 



where (> f/ > b) is used to denote samples that contain objects 
brighter than the faint or bright sample. 

In addition to the HOD parameters described above, param- 
eters to describe the spatial distribution of objects are required. 
As mentioned previously, the spatial distribution of galaxies is 
parameterized by /, a scaling factor relating the concentration 
of a DM halo to the concentration for the distribution of galax- 
ies in that halo, such that Cgai - /cdm. The three clustering 
measures - the faint autocorrelation, the bright autocorrelation, 
and the cross-correlation - require a total of two such param- 
eters. The first, /, parameterizes the spatial distribution of the 
autocorrelated samples - i.e., the distribution of any-sized satel- 
lite in any mass halo. This comprises a large catalog of objects 
and smears out environmental differences, e.g., assuming that 
/blight - /faint- The sccond, /cross parameterizes only the distribu- 
tion of faint satellite galaxies in halos with bright central galax- 
ies, as defined in the samples. This comprises a small number 
of objects and ensures that the parameter in question measures 
only the spatial distribution of satellite galaxies in galactic halos 
with ~ L* galaxies. The narrowness of the category for /coss, 
however, means that the parameter is constrained by fewer data 
points than / and depends only on the one-halo central-satellite 
term of the cross-correlation function, a term which is only dom- 
inant in the correlation function for separations less than ~0. 1 

Mpc. 

I simultaneously fit the cross-correlation function of the 
bright and faint samples, the autocorrelation function of the 
bright sample, and autocorrelation of the faint sample, which re- 
quires eight free HOD parameters: 3 parameters - Mi, a, and 
Mcut - for the faint and bright sample each and 2 parameters, / 
and /ross, to describe the spatial distribution of objects in halos. 



5. Parameter Fitting with IVIonte Carlo IVIarkov 
Chains 

5.1. MCMC 

To estimate the best-fit parameters and constrain the errors 
in their values, I use a Monte Carlo Markov Chain (MCMC) 
method. The algorithm explores the probability distributions of 
a large number of parameters such that the distribution of points 
in parameter space visited by an infinitely long MCMC chain 
follows the underlying probability density exactly. MCMC both 
measures the best-fit parameters and constrains the errors in the 
parameter estimation and a sufficiently long chain will not be 
hampered by local minima in the parameter space. 

A MCMC chain is created by choosing points in the parame- 
ter space based on a trial distribution around the previous point in 
the chain, stepping along the eigenvectors of the covariance ma- 
trix of the parameters. Trial steps that result in a smaller are 
accepted while steps that do not are accepted with some prob- 
ability, P - e^^^ For a finite chain, when the distribution 
of steps reflects the underlying distribution with sufficient accu- 
racy, the chain is said to have converged. In addition, if a chain 
is started far outside the region of high probability, an initial sec- 
tion will be unrepresentative and must be discarded; this initial 
section is referred to as 'burn in.' This truncation is unnecessary 
or can be minimized if the chain is started from a point already 
known to lie in the region of high probability. For my chains, I 
start the MCMC at the point in parameter space that has been 
identified as the likelihood maximum using a simple downhill 
simplex minimization algorithm for the data. For good mea- 
sure, the trial steps are examined for large-scale modes indicative 
of burn-in and convergence is tested by calculation of the power 
spectrum of the chain as is described bV iDunklev et al., (,2005i) . 

5.2. Testing Parameter Estimates witli Mock Gaiaxy Cataiogs 

To test the constraining power of my sample using this method, 
I create test catalogs, populating halos from a numerical simu- 
lation with reasonably-chosen values for the HOD parametersQ 
Analyzing this mock data and comparing the derived HOD pa- 
rameters to the input will show how well the analysis can work 
and how good the data needs to be for the analysis to provide 
robust constraints. The HOD parameters are listed in Table |2] 
For the spatial distribution, both f - I and /ross = 1 ■ MCMC is 
used to find the best-fit values of the eight free parameters and 
characterize the error distributions. 



Table 2. HOD Parameters for Mock Data Sets 



Catalog 






Ml 


a 




bright 
faint 


1.4 X 10" 

2.5 X 10" 


9.9 X 10'" 
4.9 X 10" 


2.0 X 10'" 
6.7 X 10'- 


1.0 
0.9 


2.0 X 10'" 
3.5 X 10" 



Note - All masses are in h Mq. 



The simulation use d here was pe r forme d using the Hashed 
Oct- Tree code of Warren & Salmonl (Il993b, and is similar to 
those presented in Seljak & Warren (2003) and IWarren et al.l 
(12006'). The simulation box is 400/i ' Mpc on a side with 1280^ 
particles and a particle mass of 2.5 x lO'^/i ^'MfT). Halos are iden- 
tified by the friends-of-friend technique dDavis et al.lll985h with 



* In this case, values are chosen to correspond roughly to the best-fit 
parameter values found for the observed sample. 
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a linking length of 0.2 times the mean interparticle separation. 
For galaxies in my faint sample, the corresponding halo in the 
simulation has roughly ~100 particles. The simulation uses a 
cosmology of Q,„ + Q.^ = 1, O,,, = 0.3, erg - 0.9, h - 0.7, and 
Qi = 0.04; I use the output at the redshift of z = 0.0. 

I create two volumes of mock data. In the first, I use the full 
simulation box to look for biases in finding the best-fit values for 
the HOD parameters and to describe the degeneracies between 
parameters. In the second, I use a volume of data similar to the 
volume of available data to show how well the analysis will work 
on the real data and how robust the the constraints will be. I also 
use this volume of data to test systematics in the SDSS data and 
the ability of the method to recover different values for / and 

/cross ■ 

5.3. Tests Using Idealized Mock Data Set 

The full simulation box of 400 Mpc on a side is 42 times 
the volume of the volume-limited SDSS data set. For this box, I 
calculate the bright and faint galaxy autocorrelations, the cross- 
correlation, and all associated covariance matrices using 400 
jackknife samples. Using the full simulation box, the diagonal 
eiTors on the faint autocorrelation and the cross-coiTelation func- 
tions are significant ly less than 5%. T he errors on the dark matter 
power spectrum (Smith et al. 2003), however, are likely larger 
than 5%. Analyzing the full simulation box, then, will only con- 
strain the errors on the dark matter power spectrum and say noth- 
ing about the biases in and degeneracies between HOD parame- 
ters. In order to evade this problem, I replace the measured faint 
autocorrelation and the cross-coiTelation with the exact values 
for the chosen HOD. 



Table 3. Unmarginalized, Best-Fit Parameters for Mock Data 
Sets 



Parameter 


Input Value 


Best-Fit with 68% errors 
Idealized 115 /i"' Mpc box 


faint: 

log(M,) 


12.826 


12.836^^ 


12 673+" "^ 


or 


0.9 


0.900!-- 


n gt^+0.124 
-^-0.094 


log(Mcu,) 


11.544 


11 561+""" 


11.691!»ij^ 


bright: 

log(M,) 


13.301 


n ^4.1 +0"'*' 


13.304!;j;^3i 


a 


1.0 


1 042+"""" 


1 055+"'-^ 


log(Me„t) 


12.301 


12 157+''-223 


10.888!l;«^ 


/ 


1.0 


1 017+"'" 
^•"^ -0.159 


0.890!-- 


/cross 


1.0 


0.981!SJ;ii 


1 051+'"**^ 


Note - Masses are in h 'Mq. 



imations in the calculation of the coiTelation function, such as 
use of spherical halo exclusion in place of elliptical exclusion. 

Degeneracies between pairs of parameters are investigated in 
Fig.|2] In each panel, the 68% and 90% confidence intervals are 
shown, marginalizing over all other parameters. / and /cross are 
more degenerate with the faint HOD parameters than the bright 
parameters. The faint sample has many more objects than the 
bright sample, so its constraining ability is greater In addition, 
/cross is only constrained by the one-halo, central-satellite term 
of the cross-correlation, which depends only on the faint HOD 
parameters. However, the degeneracies between HOD parame- 
ters - Ml, a, Mcut - in a single sample are much larger Finally, 
the degeneracy between / and /cross is pronounced and behaves 
such that larger values of / are accompanied by larger values of 
/cross and vice versa. 

The behavior of these degeneracies can be seen from the 
form of the HOD itself. Compared to an arbitrary satellite func- 
tion, {Nsai), increasing the parameter Mi reduces the overall 
number of satellites. This, however, can be compensated at 
masses larger than Mi by increasing the slope of the function 
a. At smaller masses, increasing Mi and a reduces the fraction 
of central-satellite pairs. Decreasing the roll-off parameter, Mcut, 
marginally increases the fraction of central-satellite pairs. Thus, 
larger Mi indicates larger a and smaller Mcut. In relation to the 
parameter /, increasing a tilts the correlation function to pairs in 
halos with larger masses and smaller concentrations; increasing 
/ and /cross counteracts this by boosting the concentrations. 

5.4. Tests Using SDSS-Sized Mock Data Set 

A subsample of the simulation box that is 1 15 /z^' Mpc on a side 
is equal in volume to the SDSS data sample. Using 196 jackknife 
samples (similar to the 205 samples used in the data) and esti- 
mating all correlations and covariances from the simulation, the 
size of the error bars are larger than in the case of the idealized 
box. 

The results of the parameter estimation are shown in Table 
[3] and are particularly striking in the case of the roll-off param- 
eter, Mcut, where the parameter space of likely values extends 
to Mcut < Mmin. Once Mcut drops below the minimum mass for 
central galaxies, it no longer has any effect on the coiTelation 
function, resulting in extremely large eiTors on the parameter es- 
timation. At the volume of the SDSS data set, f - I can be 
distinguished from / = 2, but values of /cross - 0.5 - 2 are con- 
sistent. 

The resulting eiTors in the parameter estimation are not di- 
rectly applicable to the SDSS data. In addition to the issue 
of sample variance, the populated simulations are an idealized 
dataset where the halos are exactly populated according to a 
HOD, there is no variance in the concentration of halos, and the 
cosmology of the populated simulation is the same as that used 
in the model. The presented analysis, however, demonstrates ap- 
proximately the degree to which the HOD and the radial distri- 
bution of satellites can be constrained with the current volume- 
limited SDSS sample. 



In Table [3] the derived best-fit parameters are listed along 
with their unmarginalized 68% eiTors, constraining / and /cross 
to within 10-15%. The table shows that the true values of all the 
parameters fall within the 68% confidence intervals. Deviations 
from true parameter values may be attributed in part to degen- 
eracies between parameters, to sample variance, or to approx- 



5.5. Tests Including Systematics and Varying f and /cross 

Using the SDSS-sized mock data set, I include the effects of the 
systematics in the SDSS data: 1.) issues with the sky subtrac- 
tion (see Section O and 2.) incompleteness in the spectroscopic 
catalog (see Section|2|. 
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Fig. 2. The 68 and 90% marginalized confidence intervals are shown for pairs of parameters in the idealized data set. All masses are 
shown in log(h^^ Mq). Crosses indicate the input parameter values. 



As discussed previously, the sky subtraction problem is angle 
dependent, over-predicting pairs at separations less than 20-30" 
and under-predicting by ~ 5% pairs at separations between 
40 - 90". At the median redshift of the data, 30" is 16/! ' kpc. 
The overestimate in pairs at small angular separations is uncon- 
strained, but in the data ~ 25% of the objects in the first bin of 
the correlation function (from 10 to 16/i ' kpc) are found at sep- 
arations less than 20". I assume that all of the 25% are spurious 
pairs and add that fraction of pairs into the mock data at sepa- 
rations less than 16/;"' kpc. 40 - 90" is 22 - 50/;"' kpc at the 
median redshift of the data. In the mock data, I remove 5% of 
the pairs with separations between 22 and 50/;"' kpc. 

In addition to including these refinements, I test different 
values of / and /cross in order to show that a range of values 
can be reproduced. The catalogs produced are MAI (/ = 0.5, 
/cross = 1); MBl (/ = 1, /cro,s = 1.5); and MCl\f = 1, 
/cross = 0.5). I also create a sample MA2 where the global com- 
pleteness is exactly 90%. This overstates the effect of incom- 
pleteness in the data where I use all areas with completeness 
better than 90%. 

Results are shown in Table|4] While the errors in these cases 
can be larger than the eiTors without accounting for systematic 
effects, there are no biases in the estimated HOD values and a 
full range of values for / and /cross are recoverable from the data. 



6. Results 

The analysis applied to the test samples can be repeated for the 
SDSS samples. First, the values of Mn,ax(cen) for the bright and 
faint sample must be estimated. With those values, the fiducial 



Table 4. Unmarginalized, Best-Fit / and /cross for Mock Data 
Sets 



Mock Sample / /c 





Input 


Best-Fit 


Input 


Best-Fit 


MAI 
MBl 
MCI 
MA2 


0.5 

1 

1 

0.5 


03,-7+0.314 

"■-"^ ' -0.179 

0.852!«;«^ 
0.987!«;^^^ 


1 

1.5 
0.5 
1 


1 021+'"'''' 

n s 1 a +0.960 
1 091+'-^" 



Note - Errors are 68% confidence intervals. 



results for the HOD parameters can be found. Then, the robust- 
ness of the best-fit values of / and /cross is tested. Finally, the 
results are compared to weak lensing results. 

6.1. Estimating M^^^^^en) 

The maximum halo mass, Mmax(cen), for the bright and faint sam- 
ples must be measured from additional SDSS samples as de- 
scribed in Section HI A luminosity-threshold sample of M, < 
-21 is used to calculate MmaxCcen) for the bright sample and a 
luminosity-threshold sample of M,- < -19 is used for the faint 
sample. For the Mr < -19 sample, I fit 4 parameters. Mi, a, 
Mcut, and /. For the Mr < -21 sample, the number of objects 
with projected separations less than = O.l/i"' Mpc is small. 
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SO constraints on / are extremely difficult to obtain, and I set 
f - 1 and fit only 3 parameters. 



Table 5. HOD Parameters for SDSS Luminosity-Threshold 
Samples 



Parameter 


Mr < - 


-21 


M,- < - 


-19 




-^max{cen) 


1.03 X 


10" 


5.44 X 


10" 


1 


M, 


1.93 X 


10'" 


6.09 X 


10'2 




a 


1.25 




0.97 








7.24 X 


10^ 


1.78 X 


10" 


a. 


/ 




10'2 


0.40 


10" 


U 


68% MnaxCcen) minimum 


9.45 X 


3.92 X 


a. 


68% MmaxCcen) maximum 


1.13 X 


10'3 


6.33 X 


10" 





Note - All masses are in h Mq. 



The best-fit values and the smallest and largest values for 
Mtnax(cen) for the unmarginalizcd parameter space that contains 
68% of the values estimated by MCMC are listed in Table |5] 
Even using luminosity-threshold samples, however, the range 
of values for the bright and faint Mmax(cen) is rather large, with 
some ~20-30% variation in the faint sample and less than ~ 10% 
variation in the bright sample. Previously published results us- 
ing different volumes of S DSS data find simila r, but not iden- 
tical results. For example, IZehavi et alj (l2005l) - which fixes 
/ = 1 - find for their Mr < -21 sample that Mmaxccen) - 
5.25 X 10i2/2-iMo, Ml = 1.23 X IO^V^Mq, ando- = 1.39. For 
their Mr < -19 sample, they find Mn,ax(cen) - 3.89x lO^'/i^'Mg, 
Ml = 8.7 X 10' VMq, and a = 1.08. 

6.2. Fiducial Results for f and /c,oss 

Using the best-fit values for Mmax(cen), the best-fit HOD param- 
eters for the luminosity-binned samples are found and shown 
in Figure [3] along with the best-fit cross-correlation. The un- 
marginalized best-fit parameters are listed in Table |6] The 



Table 6. Unmarginalized, Best-Fit Parameters for SDSS Data 

Using Best-Fit Mmax(cen) 



Parameter 


Best-Fit with 68% errors 
all bins omitting first bin 


faint: 

log(M,) 


12 75+°-"' 

'-^-0.09 


12 75+" ''' 

^^•'-'-0.10 


or 




90+° °^ 

"•'"-0.06 


log(M,„t) 


9 09+^ '* 

'■"'-0.09 


9 48+'-^* 

'•^°-0.48 


bright: 

log(M,) 


13 29+° '- 


13 29+" '^ 

^■'•^'-0.26 


a 


96+" '^ 

"•'"-0.13 


99+'' '2 

"•"-0.26 


log(M,„t) 


10 21+2" 


10.996!;|^ 


/ 


41+"" 

"■^^-0.17 




ycross 


"■'•^-0.41 


0.83tiB 


Note - Masses are in h 'Mq. 



^ 1000 



100 



10 =- 



0.01 




Tp (h"' Mpc) 



100 



10 =- 




0.1 - 



1012 1013 

M (h-' 



1014 1016 



Fig. 3. Top: The estimated cross-correlation (circles) compared 
to the halo model cross-correlation using the best-fit HOD pa- 
rameters (solid line). Bottom: The shape of the HOD for the 
bright sample (dotted) and the faint sample (solid). 



marginalized likelihoods for the parameters are shown in Figure 

El 

As expected from the tests of the mock data sets, the value 
of /cross is not well-constrained by the data, although the best- 
fit values suggest that the galaxy distribution is consistent with 
the expected DM distribution, /cross = 1 • There is a wide range 
of values from ~ .5 to ~ 2 within the unmarginalized 68% er- 
ror ellipse. Extremely small values, however, are not supported 
(/cross < 0.5). /, as expected, is better constrained, with values 
from ~ .2 and ~ 1 within the unmarginalized 68% error ellipse. 
/ < 1 is very probable. This may suggest that the distribution 
of galaxies is shallower than the dark matter distribution, which 
would still b e consistent with previous work at galaxy an d cluster 
scales (e.g.. iNagai & Kravtsovll2005l: IChen et al.ll2006l) . Finally, 
recalling the degeneracy between / and /cross, if one were to 
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BRIGHT 




Fig. 4. The normalized, marginalized likelihoods for all parame- 
ters (Ml, a, and Mem for the bright and faint samples and / and 
/cross) in the SDSS data. All masses are shown in log{h^^MQ). 



guess that the best-fit estimate of / was biased to smaller val- 
ues, one would also have to assume that the most likely value of 
/cross is also too small. 

6.3. Tests of the Robustness of the Fiducial Results 

In addition to the fiducial data set, I test a shortened data set, 
omitting the first bin in the correlation functions (physical sep- 
arations smaller than 16 /i ' kpc), with unmarginalized results 
also shown in Table|6] As discussed in Section[3] this bin is most 
affected by sky subtraction issues in the photometry. Omitting it 
does not appreciably change any of the results, aside from inflat- 
ing the error bar on /cross, the parameter most sensitive to small 
projected separations. 

Because Mmax(cen) is estimated from different luminosity 
samples than the faint and bright samples, I test the effects of the 
variances in Mmax(cen) estimates. I run the MCMC for four sam- 
ples, changing in turn the bright and faint Mn,ax(cen) to the lower 
and the upper value in the 68% unmarginalized confidence in- 
tervals and leaving the other Minax(cen) fixed to the fiducial value. 
The test samples are described and results are presented in Table 
[T] Here, the parameter estimates look very similar to the fiducial 
values, except in the case where a lower estimate of M,nax(cen) 
in the bright sample was used. In all cases, however, the errors 
overlap and all show consistent results of /ross- 

In Figure |4] the parameter Mcut has an oddly shaped and 
broad likelihood distribution (discussed briefly in the previous 
section). The marginalized likelihoods for Mcut in the SDSS 
data suggest rather small values. As noted previously, when 
Mcut < Mnin, A^ciit Contributes nothing to the HOD. One may 
be tempted to suggest that the small Mcut values are due to poor 
constraints of the data. However, small values may be well- 
motivated, at least for the bright sample. As the bright Mcut 
increases, the amplitude at small scales in the autocorrelation, 
where the one-halo central-satellite term dominates, decreases. 
The amplitude in the cross-correlation, however, is unaffected, 
since the central-satellite term in the cross-correlation samples 
the part of the faint sample satellite distribution that is unaffected 
by Mcut. Thus, large values of Mcut for the bright sample would 
cause the bright autocorrelation and the cross-correlation to in- 
tersect, which is not supported by the data (see Fig.[Tli. 



Table 7. Unmarginalized, Best-Fit Parameters for SDSS Data 
Varying M^ax (cen) 



faint: 


log(Mi) 


1 '7C+0.16 

iZ. / J_(J Qg 


1 o nc+O 19 

12.95_o,5 


1 n T/C+0 n 

12-76_oio 


1 o ac+O.lS 

12.75_oi^ 




or 




1 00+"'^ 

-0 10 


91+0.06 

-0 06 










9 01+2 "^ 


q 04.+ 1-34 


q ^4+1 69 
^■-'^-0.34 


bright: 

log(M,) 


13 29+" '^^ 

-0.18 


13 52+''-''* 




13 31+" " 




a 


"■^'-0.14 


1 17+"" 




1 00+" '" 

^^■""-0.14 




log(M,„,) 


11.34!V«5 


9 70+^ 

^- '"-0.70 




Q Q4+2 51 
^■^^-0.94 


/ 






92+" **' 


43+""' 

"•^-^-0.20 


0.47+JJ;" 


Across 




89+' 


QQQ+1.87 


1 00+ '-2** 

^•""-"-0.49 


0.88:i;» 



Note - All masses are in A^'Mq. 

" For the bright sample, M„,„(cc„) = 1.13 x IO'^/z^'Mq. 

For the bright sample, Mmax(cen) = 9.45 x IO'^/j^'M©. 

For the faint sample, M„,„(ccn) = 6.33 x IO'^Ii^'Mq. 

For the faint sample, M„,ax(ccn) = 3.92 x 1O"/!"'M0. 



6.4. Comparison to Weak Lensing Results 

If satellites are found around galaxies with the same distribu- 
tion as the dark matter profile, the cross-correlation of bright 
and faint galaxies may be related to the galaxy-mass cross- 
correlation. The projected galaxy-mass cross-correlation is pro- 
portional to the weak lensing measurement of the surface density 
contrast, AZ(rp). Given the larger errors in the weak lensing mea- 
surements at the scales at which the satellite distribution and the 
dark matter distribution can be compared, a comparison of sim- 
ple power-law fits may be sufficient, assuming a pure power-law 
density profile for X(rp) at and below those scales. 

In comparison to the weak lensing measurement, the cross- 
correlation of bright and faint galaxies has a somewhat steeper 
overall p ower-law slope than t he weak-lensing measuring. For 
example, ISheldon et al.l (|2004|) find a power-law slope of y = 
-0.76 + 0.05 for projected radii between 0.025 and 10 /z"' 
Mpc and their lens sample is peaked near an r-band magni- 
tude of -21 . In add ition, comparisons to the measurements from 
iMandelbaum et all ( f2006bl) for a lens sample with r-band mag- 
nitudes between -21 and -20 show similar power-law slopes, 
y = -0.79+0. 04 forprojectedradiibetweenO. 02 and2/!"' Mpc. 
Both measurements are shallower than the power-law slope of 
the cross-correlation of bright and faint galaxies, which show 
y - -0.86 + 0.02 at projected radii between 0.01 and 16 h^^ 
Mpc0 

The cross-correlation function, however, is more sensitive to 
the distribution of satellites around central galaxies at smaller 
separations. The best-fit power-law slope for the projected cross- 
correlation remains fairly consistent, with y - -0.88 + 0.04 
for projected radii between 0.01 and 0.25 /i ' Mpc and y = 
-0.93 + 0.07 for projected radii between 0.01 and 0.1/j ' Mpc. 
IMandelbaum et all ( l2006bl) shows larger error bars at small radii. 



' The power-law slope y for the cross-correlation is calculated by 
fitting for Wp(rp) = Ax rj,. 



Sample A° B * C^ D '' 
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with rough consistency with my results, y = -0.98+O.lOforpro- 
jected radii between 0.02 and 0.2 /i"' Mpc and y = -0.77 + 0.20 
for projected radii between 0.02 and 0.1 /; 'Mpc. The best-fit 
power-laws for the cross-correlation of bright and faint galaxies 
and the galaxy-mass cross-correlation around bright galaxies are 
statistically consistent at the most significant separations. This is 
consistent with a result of /d-oss = 1 - the satellite distribution 
follows the DM halo distribution. 

7. Conclusions 

The spatial distribution of satellite galaxies around ~ L, galax- 
ies can illuminate the relationship between satellite galaxies 
and dark matter subhalos and aid in developing and testing 
galaxy formation models. The projected cross-correlation of 
bright galaxies with faint galaxies offers a promising avenue to 
putting constraints on the radial distribution of satellite galaxies. 

Results from SDSS data of 5104 deg^ on the sky suggest 
that the spatial distribution of galaxies, parameterized by / and 
/cross, is generally consistent with the dark matter distribution 
for halos in a flat ACDM cosmology with = 0.3, Ob - 0.04, 
o"g = 0.9, and h = 0.7. /cross parameterizes the distribution of 
'faint' satellites galaxies around 'bright' central galaxies. There 
are no strong constraints on the value of /ross, however, with a 
wide range of values from ~ .5 to ~ 2 in the unmarginalized 
68% error ellipse, where /cross - 1 would indicate a distribution 
equivalent to the DM halo profile. /, on the other hand, parame- 
terizes the distribution of both bright and faint satellite galaxies 
around central galaxies of all sizes. Constraints on the value of / 
are significantly better - ruling out f -2 - and are suggestive of 
somewhat smaller concentrations for galaxy distributions com- 
pared to the concentration of dark matter halos (/ < 1). The spa- 
tial distribution of galactic satellites remains poorly constrained 
and future studies in this field will require a greater volume of 
data. 

The results presented here use a flat ACDM cosmol- 
ogy with cTg - 0.9 which is consistent with the results of 
the W ilkinson Microwave Anisotropy Probe (WMAP) 1-Year 
Data dSpergel et all |2003|) . Recent results fror n the WMAP 
5-Yea r Data supports a smaller erg - 0.8 jDunklev et alj 
120081) . iDuffv et all (l2008h claim that halo concentrations using 
a WMAP5 cosmology are lower than in a WMAPl cosmology 
by 16-23%. If halos are less concentrated than the model pre- 
dicts, the fit can compensate by decreasing /. Since / is defined 
relative to the profile of the DM halo, this implies that the true 
/ and /cross are actually larger than the estimated / and /cross- 
The difference in halo concentration is, however, fairly modest 
compared to the eiTor estimate in / and /ross- 

The results of this study have implications for previous 
works on constraining the satellite spatial distribution, galaxy 
clustering models, and weak lensing studies of the galaxy-mass 
cross-co rrelation. For exam ple, the results are consistent with 
those of IChen et al.l (|2006): satellite galaxies are distributed in 
a way that is consistent with or l ess concentrated th an the dark 
matter distribution. The results of IChen et al.l ( l2006l) apply only 
to isolated galaxies, while the cross-correlation approach applies 
to any parent galaxy. This suggests that the environment of par- 
ent galaxies is relatively unimportant. My result is also gen- 
erally consistent with cluster-sized simulations which include 
cooling and star formation. In such simulations, while the dis- 
tribution of dark matter subhalos is less concentrated than the 
dark matter, the stellar components of satellit e galaxies more 
close l y follow the dark matter distribution ( Nagai & Kravtsovl 
120051; iMaccio et ani2006l) . However, the possible discrepancy in 



values for / and /cross suggest that the spatial distribution of a 
class of satellite objects may depend on the size of the host halo: 
gal axy-sized ha los, groups and clusters. 

IChenI (I2OO8I) suggest that the projected radial distribution of 
satellite galaxies depends upon satellite color, such that redder 
satellites have a significantly steeper density profile than bluer 
satellites. A cross-correlation analysis with a different selection 
function for satellites then should show such a color dependence. 
In order to achieve better constraints on the spatial distribution 
of satellite galaxies and its environmental and color dependency, 
larger data sets are required. One possibility for increasing the 
volume of data is to apply a similar analysis to the photometric 
redshift catalog, instead of limiting the sample to objects with 
spectroscopic redshifts. The SDSS photometric redshift catalog 
goes far deeper than the spectroscopic catalog and without the 
spectroscopic catalog's fiber collision problem. However, pho- 
tometric redshifts are less accurate and may complicate the sta- 
tistical eiTor analysis. Given that future large surveys such as 
the Dark Energy Survey and the Panoramic Survey Telescope & 
Rapid Response System (Pan-STARRS) which will produce sig- 
nificantly more photometric redshifts than currently available, it 
would be interesting to perform a cross-correlation analysis with 
photometric redshifts. 
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